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Abstract. Quantized chaotic systems are generically characterized by two 
energy scales: the mean level spacing A, and the bandwidth oc h. This implies 
that with respect to driving such systems have an adiabatic, a perturbative, 
and a non-perturbative regimes. A "strong" non-linearity in the response, due 
to a quantal non-perturbative effect, is found for disordered systems that are 
described by random matrix theory models. Is there a similar effect for quantized 
chaotic systems? Theoretical arguments cannot exclude the existence of an 
analogous "weak" version of the above mentioned non-linear response effect, but 
our numerics demonstrates an unexpected degree of semiclassical correspondence. 



1. Introduction 

1.1. The two energy scales in Quantum chaos 

The name "Quantum Mechanics" is associated with the idea that the energy is 
quantized. For generic (chaotic) system the mean level spacing is A oc fi d , where 
d is the dimensionality of the system. However, one should recognize that there is 
a second energy scale Af, oc h which is introduced by Quantum Mechanics. This h 
energy scale is related to the chaos implied decay of the classical correlations. It is 
known in the literature as the "non- universal" energy scale [I] , or as the "bandwidth" 
|2] • The dimensionless bandwidth is defined as b — A^/A. For reasonably small Ti one 
has b > 1. 

This observation, of having two energy scales, has motivated the study of Wigner 
model [3] within the framework of random matrix theory (RMT). This model, which 
is defined in terms of A and Ah, is totally artificial: it does not possess any classical 
limit. Still note that it can be re-interpreted as a model for the motion of a particle 
in a quasi one-dimensional disordered lattice 0] . 

The main focus of Quantum chaos studies (so far) was on issues of spectral 
statistics [5]. In that context it turns out that the sub-?i statistical features of the 
energy spectrum are "universal" , and obey the predictions of random matrix theory. 
Non universal (system specific) features are reflected only in the large scale properties 
of the spectrum (analyzing energy intervals > A&). 

1 . 2. Regimes in the theory of driven systems 

In recent years we have made some progress in understanding the theory of driven 
quantized chaotic systems jH] Driven systems are described by Hamiltonian 
H(Q, P., x(t)), where x(t) is a time dependent parameter. Such systems are of interest 
in mesoscopic physics (quantum dots), as well as in nuclear, atomic and molecular 
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physics. The time dependent parameter x(t) may have the significance of external 
electric field or magnetic flux or gate voltage. Linear driving x(t) = Vt is characterized 
by one parameter (V), while more generally a periodic driving x(t) = Af(t) is 
characterized by both amplitude (A) and frequency (CI). Due to the time dependence 
of x(t), the energy of the system is not a constant of motion. Rather the system makes 
"transitions" between energy levels, and therefore absorbs energy. 

The main object of our studies is the energy spreading kernel Pt(n\m). Regarded 
as a function of the level index n, it gives the energy distribution after time t, 
where m is the initial level. Having two quantal energy scales (A, A&) implies the 
existence of different quantum-mechanical (QM) V regimes jS], or more generally 
(A, fi) regimes [7], in the theory of Pt(n\m). Most familiar is the QM adiabatic regime 
(very very small V), whose existence is associated with having finite A. Outside of the 
adiabatic regime we are used to the idea that there is a perturbative regime, where the 
Fermi golden rule applies, leading to a Markovian picture of the dynamics, with well 
defined transition rates between levels. Less familiar |SJ[7] is the QM non-perturbative 
regime (V is quantum mechanically large, but still classically small) whose existence 
is associated with the energy scale A&. As implied by the terminology, in the QM 
non-perturbative regime perturbation theory (to any order) is not a valid tool for the 
analysis of the energy spreading. Consequently the Fermi golden rule picture of the 
dynamics does not apply there. 

1.3. Linear response theory 

Of special importance (see discussion below) is the theory for the variance 8E(t) 2 = 
^2 n Pt{n\m)(E n — E m ) 2 of the energy spreading. Having 6E(t) oc A means linear 
response. If SE(t)/A depends on A we call it "non-linear response". In this paragraph 
we explain that linear response theory (LRT) is based on the "LRT formula" 

/OO 1 
=-F t {u)C(u) (1) 

Two spectral functions are involved: One is the power spectrum C(u>) of the 
fluctuations, and the other F t (Lo) is the spectral content of the driving. See Eq.(£lJ 
and Eq.JSJ) for exact definitions. A special case of Eq.Q) is the sudden limit (V = oo) 
for which f(t) is a step function, hence F t (uj) = 1, and accordingly 

SE = \/C(0) x A ["sudden" case] (2) 

Another special case is the response for persistent (either linear or periodic) driving. 
In such case the long time limit of F t {ui) is linear in time [e.g. for linear driving 
(f(t) = t) we get F t (ui) = t x 2wS(u)]. This implies diffusive behavior: 

8E{t) = y/2D E t ["Kubo" case] (3) 

In the latter case the expression for De oc A 2 is known as Kubo (or Kubo-Greenwood) 
formula, leading to a fluctuation-dissipation relation [S]. 

The LRT formula Eq. has a simple classical derivation 6 . From now on it 
goes without saying that we assume that the classical conditions on (A, Q) for the 
validity of Eq.QJ are satisfied (no ft involved in such conditions). The question is 
what happens to the validity of LRT once we "quantize" the system. Can we trust 
Eq.QJ for any (A, O)? Or maybe we can trust it only in a restricted regime? In 
previous publications IE], we were able to argue the following: 
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(A) The LRT formula can be trusted in the perturbative regime, with the exclusion 
of the adiabatic regime. 

(B) In the sudden limit the LRT formula can be trusted also in the non-perturbative 
regime. 

(C) In general the LRT formula cannot be trusted in the non-perturbative regime. 

(D) The LRT formula can be trusted deep in the non-perturbative regime, provided 
the system has a classical limit. 

For a system that does not have a classical limit (Wigner model) we were able to 
demonstrate that LRT fails in the non-perturbative regime. Namely, for Wigner 
model the response 5E(t)/A becomes A dependent for large A, meaning that the 
response is non-linear. Hence the statement in item (C) above has been established. 
We had argued that the observed non-linear response is the result of a quantal non- 
perturbative effect. Do we have a similar type of non-linear response in case of 
quantized chaotic systems? The statement in item (D) above seems to suggest that 
the observation of such non-linearity is not likely. Still, we argue below that this does 
not exclude the possibility of observing a "weak" non-linearity. 

2. Perturbation theory and linear response 

The immediate (naive) tendency is to regard LRT as the outcome of first order 
perturbation theory (FOPT). In fact the regimes of validity of FOPT and of LRT 
do not coincide. On the one hand we have the adiabatic regime where FOPT is valid 
as a leading order description, but not for response calculation (see further details 
below). On the other hand, the validity of Eq.QJ goes well beyond FOPT. This leads 
to the (correct) identification [7] of what we call the "perturbative regime" . The 
border of this regime [in {A, Q) space] is determined by the energy scale Ah, while A 
is not involved. Outside of the perturbative regime we cannot trust the LRT formula. 
However, as we further explain below, the fact that Eq.Q is not valid in the non- 
perturbative regime, does not imply that it fails there. 

We stress again that one should distinguish between "non-perturbative response" 
and "non- linear response". These are not synonyms. As we explain in the 
next paragraph, the adiabatic regime is "perturbative" but "non-linear" , while the 
semiclassical limit is "non-perturbative" but "linear" . 

In the adiabatic regime, FOPT implies zero probability to make a transitions to 
other levels. Therefore, to the extend that we can trust the adiabatic approximation, 
all the probability remains concentrated in the initial level. Thus, in the adiabatic 
regime, Eq.([Q) is not a valid formula: It is essential to use higher orders of perturbation 
theory, and possibly non-perturbative corrections (Landau-Zener j2j), in order to 
calculate the response. Still, FOPT provides a meaningful leading order description of 
the dynamics (i.e. having no transitions), and therefore we do not regard the adiabatic 
non-linear regime as "non-perturbative" . 

In the non-perturbative regime the evolution of Pt(n\m) cannot be extracted from 
perturbation theory: not in leading order, neither in any order. Still it does not 
necessarily imply a non-linear response. On the contrary: the semiclassical limit is 
contained in the (deep) non-perturbative regimejjj. There, the LRT formula Eq.Q is 
in fact valid. But its validity is not a consequence of perturbation theory, but rather 
the consequence of quantal- classical correspondence. 
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3. The quest for non perturbative response 

As stated above, an effect of non-linear response due to the quantum mechanical non- 
perturbative nature of the dynamics, has been demonstrated so far only for Wigner 
model [7|. There, its existence is related to the disordered RMT nature of the model 
(see discussion below). Semiclassical correspondence considerations seem to exclude 
the manifestation of this disorder-related non-linearity in case of quantized chaotic 
systems. In this Letter we explain that this does not exclude the possibility of having 
a "weak" version of this effect. We also report the results of an intense numerical 
effort aimed in finding a "weak" non-linearity in the case of a simple low-dimensional 
quantized chaotic systems. To our surprise, an unexpected degree of semiclassical 
correspondence is observed. 

It is appropriate here to clarify the notions of "weak" and "strong" effects. In 
the literature regarding the dynamics in disordered lattices one distinguishes between 
"weak" and "strong" localization effects. The former term implies that while the 
leading behavior is classical (diffusion), there are "on top" quantum mechanical 
corrections (enhanced return probability). In contrast to that the term "strong" 
implies that the classical description fails even as a leading order description. In 
the literature regarding quantum chaos we have the effect of "scarring" , which should 
be regarded as "weak" effect. "Strong" quantum mechanical effects (e.g. dynamical 
localization in ID kicked systems ^2]) are non-generic: The leading order behavior 
of generic quantized chaotic systems is typically classical. In the present context of 
driven systems, we use the terms "weak" and "strong" in the same sense: The adjective 
"weak" is associated with the (conjectured) non-linear response of quantized driven 
chaotic systems, while the adjective "strong" is associated with the (established) non- 
linear response in the corresponding RMT (Wigner) model. 

4. The numerical findings 

How do we detect non-linear response? The most straightforward way is to fix 
the pulse shape f(t) and to plot 5E/A versus A. A deviation from constant value 
means "non-linear" response. The simulations below are done for a quantized chaotic 
system. Due to obvious numerical limitations we will consider the response to one- 
pulse driving, rather than to persistent (periodic) driving. The central question is 
whether the observed non-linear effect is of semiclassical origin, or of novel quantum 
mechanical origin. We deal with this issue below. 

In our numerical simulations (Fig.l) we have considered a particle in a two 
dimensional (d = 2) anharmonic well. This model (with deformation parameter 
x = const) is defined in ^| In the energy region of interest (E ~ 3), the 

classical motion inside the two dimensional well (2DW) is chaotic, with characteristic 
correlation time r c ; ~ 1. For the following presentation it is enough to say that the 
quantum mechanical Hamiltonian is represented by a matrix TL = E + x(t)H, where E 
is a diagonal matrix with mean level spacing A w 4.3 x h d , and B is a banded matrix. 
The bandwidth (in energy units) is = 2nh/ 't c i. The bandprofile (see Fig. 2 of \W\) 
is described by a spectral function which is defined as follows: 



with implicit average over the reference state m. The bandprofile, as defined above, 
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can be determined from the the classical dynamics. This means that C(u>) » C c1 (uj) 
where C c1 (lo) is the Fourier transform of a classical correlation function C c1 (t). The 
h dependence of C(u>) is relatively weak. 

The driving pulse in our numerical simulations has a rectangular shape. This 
means /(0) = f(T) = and /(0 < t < T) = A, where T = 0.375. The spectral 
content of the driving is defined as: 

2 



F t (w) = 



f{t')e iut dt' 



(5) 



The spectral content of the driving after a rectangular pulse is F t (oj) — |1 — e lulT \ 2 . 
We have also made simulations (not presented) with a driving scheme that involves a 
positive pulse +A followed by a negative pulse —A, with the intention of considering 
eventually a persistent (multi cycle) periodic driving. However, we have realized that 
all the relevant physics is observed already in the single pulse case. Note that the 
regime diagram for either linear or (as in the following simulations) rectangular driving 
pulse, is greatly simplified, because the driving is characterized by only one parameter 
(V in the former case, A in the latter case). 

Let us look carefully at the results of the 2DW simulations (Fig.l). For small A 
we see as expected "linear response" meaning SE/A — const, as implied by Eq.QJ. 
Note that the "constant" has a weak fi dependence (a 10% effect). This is due to the 
above mentioned weak dependence of C(oj) on Ti. So this quantum-mechanical effect is 
quite trivial, and has a simple explanation within LRT. Now let us look what happens 
for large A. We clearly see a 2% deviation from linear response. In what follows we 
discuss the reason for this non-linear effect. 

For sake of comparison we also perform simulations with an effective RMT model 
that corresponds to the 2DW model Hamiltonian. The effective RMT model is 
obtained by randomizing the signs of the off-diagonal elements of the B matrix. The 
effective RMT Hamiltonian has the same bandprofile C(uj) as the original (2DW) 
Hamiltonian. Therefore, as far as LRT Eq.Q) is concerned, the response should be 
the same. Still we see that at the same A regime, as in the case of the 2DW simulations, 
we have deviation from linear response. However, this non-linear deviation is much 
much stronger. 

Looking at the curves of Fig.l, it is very tempting to regard the observed non- 
linear 2% effect in the 2DW simulations as a "weak" version of the "strong" effect 
which is observed in the corresponding RMT simulations. However, the careful 
analysis below indicates that apparently this is not the case. 



5. Discussion and analaysis 

In analyzing the validity of the LRT formula, it is instructive to consider first the 
sudden limit Eq.Q. This limit has been studied in JHj- The spreading profile P(n\m), 
after the sudden change in x, depends on the amplitude A of the perturbation. [We 
omit the time index t, which is of no relevance in this limit] . The perturbative regime is 
A < A plt , where A plt = 2nh/ (t c1 ^JC{0)). For the 2DW simulations A plt = 5.3 x fi. In 
the perturbative regime P(n\m) has a core-tail structure (the generalization of Wigner 
Lorentzian), and the variance SE 2 is determined by the first order tail component 
of the energy distribution. For A > A prt the spreading profile P(n\m) becomes 
non-perturbative. This means that the perturbative tail (if survives) is no longer 
the predominant feature. Thus the variance is determined by the non-perturbative 
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component (the "core" ) of the energy distribution. The remarkable fact is that the 
crossover from the perturbative A regime to the non-perturbative A regime is not 
reflected in the variance (see Fig. 5 of ^UJ ) . The agreement with Eq.© is perfect. 
Taking into account the "dramatic" differences in the appearance of P(n\m), this 
looks quite surprising. In fact (see Sec. 12 of JUj) there is a simple proof ^31 that 
Eq. J2J remains exact beyond any order of perturbation theory, which means that it is 
exact even in the non-perturbative regime where perturbation theory is not applicable. 

We turn back to our simulations, where we have a rectangular pulse (rather than 
step function). Here the sudden limit does not apply, and the dynamics within the 
time interval < t < T should be taken into account. If we take an eigenstate of E 
and propagate it using E + ^4B, then we get in the classical case ballistic spreading 
followed by saturation. ["Eigenstate" in the classical case means microcanonical 
distribution]. This is true for any A. Quantum mechanically we observe in the 
2DW model simulations a similar ballistic behavior JT] , whereas in the corresponding 
RMT model there is an intermediate stage of diffusion JT] . This diffusion is of non- 
perturbative nature, and it is related to the "disorder" which is artificially introduced 
via the sign-randomization procedure. The strong non-linear response effect [J] is a 
consequence of this diffusion. 

Coming back to the 2DW model, we realize that there is no "disorder" build 
into the model, and therefore no diffusion. Still, looking at Fig.l, it is tempting to 
interpret the observed 2% non-linear deviation as a "weak" version of the strong 
non-linear effect. Moreover, regarded as such, it vanishes, as expected, in the 
deep non-perturbative regime, which had been argued on the basis of semiclassical 
correspondence considerations 0. 

In order to properly determine whether the dips in Fig.l are the result of the 
QM non-perturbative nature of the dynamics, we have rescaled the vertical axis, and 
plotted the response once (Fig. 2a) versus A, and once (Fig. 2b) versus A/%. On the 
basis of the scaling we see that the strong non-linear response in the RMT case is 
indeed the result of the quantal (^-dependent) non-perturbative effect. In contrary 
to that the ^-independent scaling in the 2DW case, indicates that the non-linear 
deviation there is of "semiclassical" rather than "quantal non-perturbative" nature. 

6. Conclusion 

Theoretical arguments cannot exclude the existence of a "weak" non-linearity in the 
response of a driven quantized chaotic system, which is due to a quantum mechanical 
non-perturbative effect. But our careful numerics, regarding a simple low-dimensional 
system, demonstrates an unexpected degree of semiclassical correspondence. Our 
findings should be regarded as the outcome of an ongoing quest, which has not ended, 
that is aimed in finding novel quantum mechanical deviations from linear response 
theory. 
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FIG.l: The response 5E/A as a result of a rectangular pulse (T = 0.375). Deviation from 
5E/A = const implies non-linear response. All the data are averaged over a number of different 
initial conditions. The simulations are done with the 2DW Hamiltonian (circles), and also with the 
associated RMT model (stars). See text for explanations. 




FIG. 2: Scaled versions of Fig.l. The vertical scaling is aimed in removing the weak fi dependence 
of the bandprofile. In (b) the horizontal scaling is aimed in checking whether the deviation from linear 
response is in fact a quantal non-perturbative effect. 



